Nuclear symmetry energy from the Fermi-energy difference in 

nuclei 

Ning Wang, 1 '* Li Ou, 1 and Min Liu 1 
1 Department of Physics, Guangxi Normal University, Guilin, 541004, P- R- China 

Abstract 

The neutron-proton Fermi energy difference and the correlation to nucleoli separation energies 
for some magic nuclei are investigated with the Skyrme energy density functionals and nuclear 
masses, with which the nuclear symmetry energy at sub-saturation densities is constrained from 54 
Skyrme parameter sets. The extracted nuclear symmetry energy at sub-saturation density of 0.11 
fm -3 is 26.2 ± 1.0 MeV with 1.5a uncertainty. By further combining the neutron-skin thickness 
of 208 Pb, ten Skyrme forces with slope parameter of 28 ^ L ^ 65 are selected for the description 
of the symmetry energy around saturation densities. 
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I. INTRODUCTION 



The nuclear symmetry energy, in particular its density dependence, has received consid- 
erable attention in recent years [1-12]. As one of the key properties of nuclear matter, the 
nuclear symmetry energy probes the isospin part of nuclear force and intimately relates to 
the structure character of drip line nuclei and super-heavy nuclei, the dynamical process 
of nuclear reactions and the behavior of neutron stars. To explore the density dependence 
of the nuclear symmetry energy from sub-saturation to super-saturation densities, various 
models and experimental observables have been proposed. On one hand, the constraints 
on the symmetry energy are investigated from heavy-ion collisions [2, 13, 14]. Some ex- 
perimental data for the isospin sensitive observables are well reproduced by using certain 
forms of density-dependent symmetry energy in microscopic dynamics calculations, such as 
in the improved quantum molecular dynamics [1, 2] and the isospin Boltzmann-Uehling- 
Uhlenbeck [4, 5] calculations. In these calculations, the temperature effect and the influence 
of the isospin-independent terms of nuclear force are self-consistently involved. It is still 
difficult to clearly obtain the information of nuclear symmetry energy at zero-temperature 
by removing the influence of the isospin-independent terms. 

On the other hand, the symmetry energy is also constrained from the properties of finite 
nuclei, such as the binding energies [15-20], the neutron skin thickness [8, 21-24] and the 
pygmy dipole resonance (PDR) [25, 26]. By analyzing the more than 2000 measured masses 
of nuclei with the help of the liquid drop formula, one can obtain the mass dependence of 
the symmetry energy coefficients of finite nuclei [15, 16, 20] and the nuclear symmetry en- 
ergy at sub-saturation densities [17, 18] based on the relation between the symmetry energy 
coefficients of finite nuclei and the symmetry energy of infinite nuclear matter [8]. The ob- 
tained symmetry energy at the saturation density and its slope parameter [17] are generally 
close to the results from PDR [26] and heavy- ion collisions [1]. In addition, the symmetry 
energy is also constrained from neutron star observations incorporates the microphysics of 
both the stellar crust and core [9, 27, 28]. A relatively smaller slope parameter of the sym- 
metry energy at the saturation density 43 < L < 52 MeV is obtained from the available 
neutron star mass and radius measurements [9]. It is known that the average density of 
a finite nucleus is smaller than the saturation density due to the surface diffuseness. The 
obtained information mainly describes the nuclear symmetry energy at the sub-saturation 
densities rather than the super-saturation densities. The symmetry energies at the satura- 
tion and super-saturation densities based on the extrapolation are still very uncertain and 
more isospin-sensitive obervations should be proposed and investigated. 

The density dependence of the symmetry energy is also extensively studied with the 
Skryme energy density functionals [29-31]. The Hartree-Fock-Bogoliubov approach with the 
Skyrme force BSkl7 [30] can reproduce the 2149 measured masses with an rms deviation 



2 



of 0.581 MeV, which is comparable to the accuracy of the new finite range droplet model 
(FRDM2012) [19]. In the FRDM2012, it is found that a slope parameter of the symmetry 
energy at normal density L ps 70 MeV can give better results with an rms deviation of 0.570 
MeV. However, one should note that the corresponding slope parameter is only L = 36 
MeV from the BSkl7 parameter set which is much smaller than the result of FRDM2012. 
It is therefore necessary to further investigate the behavior of the symmetry energy at sub- 
saturation densities from the structures of finite nuclei. 

In this work, we study the nuclear symmetry energy at sub-saturation densities from the 
Fermi energies of nuclei based on various parametrizations of the Skyrme forces. The Skyrme 
interaction, originally constructed for finite nuclei and nuclear matter at saturation density, 
is a low momentum expansion of the effective two-body NN interaction in momentum space. 
Although all Skyrme forces are usually fitted to reproduce well the saturation energy and 
density of symmetric nuclear matter, they differ significantly in other characteristics of 
symmetric and pure neutron matter, in particular their density dependence [31]. The Fermi 
energies of neutrons and protons for some doubly-magic nuclei can be measured with a 
high precision. The neutron-proton Fermi energy difference and the correlation to nucleon 
separation energies of magic nuclei are directly related to the symmetry energy of nuclei, with 
which the isospin-dependent terms of the Skyrme forces and thus the density dependence of 
the symmetry energy could be constrained. The paper is organized as follows: In Sec. II, the 
correlation between Fermi energies and separation energies of nucleons are introduced. In 
Sec. Ill, the nuclear symmetry energy at sub-saturation and saturation densities is extracted. 
Finally, a summary is given in Sec. IV. 

II. CORRELATION BETWEEN FERMI ENERGIES AND SEPARATION 

ENERGIES OF NUCLEONS 

Based on the liquid drop mass formula, the binding energy of a nucleus which is taken 
as positive value is expressed as 

Z 2 

BE (A, Z) = a v A - a s A 2 ^ - a c ^ - a sym I 2 A (1) 

by neglecting nuclear microscopic corrections. / = (N—Z)/A denotes the isospin asymmetry. 
The difference between the proton separation energy [BE(A, Z) — BE{A — 1, Z — 1)] and 
the neutron separation energy [BE(A, Z) — BE(A — 1, Z)] of a nucleus is written as, 

AS = S p - S n ~ -2a c -^ + 4a sym J. (2) 

Because the Coulomb energy coefficient a c = |^ ~ 0.71 MeV is usually well determined 
from the masses of mirror nuclei [32, 33], the value of AS is directly related to the symmetry 
energy coefficients of finite nuclei. 
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On the other hand, the single-particle energies (SPE) of a nucleus can be uniquely deter- 
mined by solving the Schrodinger equations or Hartree-Fock equations based on the single 
particle potential, under the mean-field approximation. In the Hartree-Fock theory for a 
closed-shell nucleus (A, Z) the single-particle energies for states below the Fermi surface are 
given by [34, 35] 

e p = BE*(A- 1,Z- 1) - BE(A,Z) (3) 

and 

e n = BE*(A- 1,Z) -BE(A,Z). (4) 

e will be negative for bound states. (BE* = BE — E x ) is the ground state binding energy 
minus the excitation energy of the excited states associated with the single-particle states. 
The difference between the Fermi energy of neutrons and that of protons 

Ae = e F n - < (5) 

is closely related to the proton and neutron separation energies of the nucleus. Here, the 
Fermi energy e F is defined as the energy of the highest occupied quantum state in a system 
of fermions at absolute zero temperature. 

Fig. 1 shows the calculated single particle energies of occupied states for protons and 
neutrons of 132 Sn at its ground state by using the Skyrme Hartree-Fock (SHF) model with 
the parameter set SLy7 [36]. For this neutron-rich nucleus, the Fermi energy of neutrons is 
higher than that of protons, and the calculated value of Ae is 8.1 MeV. The depth of the 
single particle potential plays a dominant role for the corresponding Fermi energy of a given 
nucleus. The difference Ae closely relates to the difference between the depth of nuclear 
potentials for neutrons and protons, 

Vn - (V p + V C ) ~ ~\~§^ + 2W- (6) 

Here, Vq — denotes the Coulomb potential of a nucleus at the central position with 

the potential radius r c » 1.3 fm. The information on the symmetry potential from the 
Fermi energy difference is of great importance for the study of nuclear symmetry energy. 

It is known that under the mean-field approximation, if the single-particle motion plays 
a dominant role for the behavior of nucleons near the Fermi surface, one expects that the 
relation AS* ~ Ae holds for the closed-shell nuclei. The experimental values of Ae for six 
doubly-magic nuclei 16 0, 40 Ca, 48 Ca, 56 Ni, 132 Sn and 208 Pb are -3.53, -7.31, 6.10, -9.47, 
8.40 and 0.64 MeV, respectively [34, 35]. The corresponding separation energy differences 
AS* for these six nuclei are —3.53, —7.31, 5.86, —9.48, 8.39, 0.64 MeV, respectively. One sees 
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FIG. 1: (Color online) Single particle potentials and the single particle energies (SPE) of bound 
states for 132 Sn with the Skyrme Hartree-Fock calculation by using SLy7 force. The dashed lines 
denote the highest SPE of the occupied states for nucleons in 132 Sn at its ground state. 

that the relation AS ~ As does hold very well as expected for these doubly-magic nuclei. To 
further test the relation, the values of Ae for 19 doubly-magic or semi-magic nuclei ( 16,22 0, 

22,42g i; 40,48,60^ 42^ 56,68,78^ 130^ 100,132,134^134^ 182 , 20 8p b) ^ systematically 

investigated by using the SHF model together with 54 commonly used Skyrme forces. The 
calculated values are shown in Fig. 2 as a function of AS. Here, the masses of unmeasured 
nuclei are predicted with the Weizsacker- Skyrme mass formula combining the radial basis 
function correction (WS3 RBF ) [37]. The WS3 RBF model can reproduce the measured 2149 
masses in AME2003 with an rms deviation of 184 keV, and the predictive power is also 
remarkable [38] (the rms deviation with respect to the 154 new masses of extremely neutron- 
rich and proton-rich nuclei listed in AME2012 [39] is only 397 keV). The red squares in 
Fig. 2 which denotes the experimental data for the six doubly-magic nucleus mentioned 
previously are quite regularly located along the green line AS* = Ae. The solid circles, open 
circles and crosses denote the results of three Skyrme forces SkSC4 [40], SLy230a [41] and 
BSkl7 [30], respectively. The error bars denote the uncertainty of the model calculations 
from 54 different Skyrme forces in which the corresponding incompressibility coefficient for 
symmetry nuclear matter is = 230±30 MeV and the saturation density po = 0.16±0.005 
fm -3 . One sees that the calculated results from the traditional 10-parameter Skyrme forces 
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FIG. 2: (Color online) Fermi energy difference as a function of separation energy difference. The 
red squares denote the experimental data for six doubly-magic nuclei [34, 35]. Others are the 
calculated results with difference Skyrme forces. 

support the relation AS ~ Ae, even for the extremely neutron-rich and proton-rich nuclei. 
The correlation between the Fermi energy difference and the separation energy difference for 
neutrons and protons could be helpful for constraining the equation of state for asymmetry 
nuclear matter. Here, we would like to emphasize that the Fermi energy difference can 
effectively remove the influence of isospin-independent terms in the nuclear forces. 



III. NUCLEAR SYMMETRY ENERGY AT SUB-SATURATION AND 



Based on the calculated Fermi energy difference Ae with the 54 different Skyrme forces, 
the average deviation 



from the m — 19 nuclei mentioned previously is calculated. Fig. 3(a) and (b) show the 
average deviation as a function of symmetry energy -E sym at the density of p = 0.09 and 0.11 
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FIG. 3: (Color online) Average deviation as a function of nuclear symmetry energy E syra at the 
density of (a) p = 0.09 fm~ 3 and (b) 0.11 fm -3 , with 54 different Skyrme forces. The solid curves 
denote the parabolic fit to the squares. The dashed line show the position of 1.19 MeV. 

fm -3 , respectively. The nuclear symmetry energy in the Skyrme energy density functional 
is expressed as 

1 

E sym (p) = - 

1 k 2 /37T 2 Y /3 2/3 1 fn 

~ ki (if) 6symp5/3 - k hi2X * + ^ (8) 

with 6 sym = 3tiXi — t 2 (4 + 5x 2 ). t , ti, t 2 , t 3 , x , x±, x 2 , x 3 and a are the Skyrme parameters. 
The squares and solid curves denote the results of 54 Skyrme forces and the parabolic fit, 
respectively. One sees that the minimal deviations are located at around £ , sym (0.09) = 
23.6 and £ , sym (0.11) = 26.2 MeV, respectively. We also note that the obtained symmetry 
energies do not change appreciately if only the six doubly-magic nuclei are involved in the 
calculation of the average deviation (a). It indicates that the Fermi energy difference is 
a useful observation for studying the symmetry energy at sub-saturation densities. The 
dashed line shows the minimal value of the parabolic curve which is 1.19 MeV. The minimal 
average deviation from the 54 Skyrme forces is about 1.0 MeV. Considering the systematic 
error of the WS3 RBF mass model which is about 0.19 MeV [37], the Skyrme forces with 
(a) < 1.19 MeV reasonably well describe the Fermi energy difference for the 19 nuclei 
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FIG. 4: (Color online) Average deviation as a function of (a) saturation density po an d (b) in- 
compressibility coefficient of symmetric nuclear matter, based on the calculations of the 54 
Skyrme forces. 

mentioned previously. 

As one of the key properties of nuclear matter, the symmetry energy is particularly 
important in modelling nuclear matter and finite nuclei because it probes the isospin part 
of the Skyrme interaction. For a sensitive observation to investigate the nuclear symmetry 
energy, the influence of isospin-independent terms in the nuclear forces should be removed 
as clean as possible. In Fig. 4, we show the average deviation as a function of nuclear 
saturation density po and the incompressibility coefficient for symmetry nuclear matter. 
One cannot evidently obtain the optimal values of po and according to the average 
deviations from the 54 Skyrme forces, since these two quantities are determined by the 
isospin-independent parts of the Skyrme interactions. We also note that the strength of the 
spin-orbit interaction Wq in the Skyrme forces does not affect the value of Ae generally due 
to the cancelation between protons and neutrons, which is helpful to remove the influence of 
the shell effect on the symmetry energy. For example, the value of Ae for 132 Sn only changes 
by 0.7% with a variation of the strength of the spin-orbit interaction by 32% according to 
the SLy7 calculations. 

From the 54 Skyrme forces, 17 parameter sets with (a) < 1.19 MeV are selected for the 
description of E sym (p). These selected Skyrme forces can well reproduce the experimental 
data for the Fermi energy difference of the six doubly-magic nuclei. We note that the 
obtained symmetry energies from these 17 Skyrme interactions are close to each other at 
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FIG. 5: (Color online) The same as Fig. 3, but at the density of p = 0.13 and 0.16 fm 3 , respectively. 

the densities around p c = 0.11 fm -3 , and the value of E sym (p c ) = 26.2 ± 1.0 MeV with 1.5a 
uncertainty. It could be much more useful if the slope parameter of the symmetry energy 
at the saturation density po 

can be also well constrained with this approach. Unfortunately, we find that the obtained 
uncertainties of the symmetry energy at densities lower and higher than p c gradually increase. 
Fig. 5 shows the average deviation as a function of nuclear symmetry energy £- S y m at the 
density of p = 0.13 and 0.16 fm -3 , respectively. The parabolic behavior of the average 
deviation becomes unclear with increasing of density and even disappears at the saturation 
density. The corresponding uncertainties of the symmetry energy at p = 0.13 and 0.14 fm -3 
increase to 1.5 and 2.0 MeV, respectively. It indicates that the slope parameter L can not 
be well constrained by using the neutron-proton Fermi energy difference uniquely, although 
the nuclear symmetry energy at sub-saturation density can be well described. 

To further analyze the density dependence of the symmetry energy at the saturation 
density, the values of the symmetry energy at p c and p , the slope parameters L, and the 
neutron-skin thickness Ai? np of 208 Pb are listed in Table I. One sees that for the seven 
Skyrme forces, i.e. SkSC4, SkSCl, v075, v080, v090, vl05, MSk3, the corresponding slope 
parameters L are very small and even negative. In Ref. [31], it is thought that the BSk, 
SkSC and MSk families under-predict both the symmetry energy and its derivative at the 
saturation density. To further check these selected Skyrme forces for the description of other 
physical quantities, we study the neutron-skin thickness of 208 Pb. The linear relationship 
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FIG. 6: (Color online) Slope parameter of nuclear symmetry energy as a function of neutron- 
skin thickness of 208 Pb. The region between the two dashed vertical lines denotes the measured 
neutron-skin thickness of 208 Pb, 0.131 < Ai? np ( 208 Pb) < 0.218 fin [23-26]. 

between the slope parameter L and the Ai? np of 208 Pb was observed in Refs.[8, 22]. Fig. 
6 shows the slope parameter L of these Skyrme forces as a function of the corresponding 
neutron-skin thickness of 208 Pb. The corresponding values of Ai? np ( 208 Pb) from the seven 
Skyrme forces with L < 10 MeV are about 0.11 fm. The recent experimental measurements 
[23-26] for the AR np of 208 Pb show 0.131 ^ Ai? np ( 208 Pb) < 0.218 fm. It implies that the 
seven Skyrme forces with small L values are not suitable for the description of the symmetry 
energy at around saturation densities. The other ten Skyrme forces with 28 ^ L ^ 65 MeV, 
i.e., Ska25s20, Ska35s20, SV-sym32, SLy2, SLy6, SLy7, SLy230a, BSkl7, Skzl and SkT7, 
reasonably well describe both the Fermi energy difference and the neutron-skin thickness of 
208 Pb. The obtained central value of L with this approach is generally consistent with the 
result from the neutron star mass and radius measurements [9]. Out of the ten forces, four 
with the smallest average deviation [(a) ^ 1.12 MeV, see Table I] have values of L = 56 ± 9 
MeV for the slope parameter, and the Skyrme forces Ska25s20, Ska35s20, SV-sym32 were 
also recommended in the recent study of Ref. [31]. 

Fig. 7 shows the calculated nuclear symmetry energy as a function of density. The squares 
(with 1.5a uncertainty as the error bar) denote the results from the 17 selected Skyrme forces 
listed in Table I. Other symbols and curvatures denote the results from the selected ten forces 
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FIG. 7: (Color online) Nuclear symmetry energy as a function of density. The squares denotes the 
results from the 17 selected Skyrme forces based on the Fermi energy difference. Other symbols 
and curves denotes the results of ten selected Skyrme forces with which both the Fermi energy 
difference and the neutron-skin thickness of 208 Pb can be reasonably well described. 

which can well describe the neutron-skin thickness of 208 Pb simultaneously. The symmetry 
energy at the saturation density from the ten selected Skyrme forces is E syra (p ) = 31. 9 ±2.1 
MeV with 1.5a uncertainty. 

IV. SUMMARY 

The correlation between the neutron-proton Fermi energy difference Ae and the separa- 
tion energy difference AS for some doubly-magic and semi-magic nuclei is analyzed with the 
Skyrme energy density functionals and nuclear masses, with which nuclear symmetry energy 
at sub-saturation densities is constrained from 54 different Skyrme forces. The experimental 
data and the Skyrme Hartreee-Fock calculations show AS* ~ Ae for closed-shell nuclei, even 
at the extremely neutron-rich and proton-rich cases. The correlation between Ae and AS 
is a good observation for studying the nuclear symmetry energy at sub-saturation densities 
which probes the isospin-dependent part of the Skyrme interaction, since the cancelation 
between protons and neutrons directly removes the influence of isospin-independent terms. 
The extracted symmetry energy from 17 selected Skryme forces at the density of 0.11 fm -3 
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TABLE I: Nuclear symmetry energy and neutron-skin thickness Ai? np of Pb with the selected 
17 Skyrme forces, (a) denote the calculated average deviation according to Eq.(7) for the 19 nuclei. 
The unit of AR np is fm, and those of others are MeV. The bold-face entries denote the four forces 
with the smallest average deviation and reasonable neutron-skin thickness of 208 Pb. 
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is about 26.2 ± 1.0 MeV. The slope parameter of symmetry energy is also studied by fur- 
ther combining the neutron-skin thickness of 208 Pb. Out of 54 Skyrme forces, ten with 
28 ^ L ^ 65 MeV can reasonably well describe both the Fermi energy difference and the 
neutron-skin thickness of 208 Pb. Within the ten forces, four Skyrme forces with the smallest 
deviation (i.e. Ska25s20, Ska35s20, SV-sym32 and SLy2) have values of L = 56 ± 9 MeV 
for the slope parameter. The structures and masses of finite nuclei at the ground states 
is helpful to obtain the information of symmetry energy at the sub-saturation densities. 
One should also note that the uncertainty of the symmetry energy extracted from nuclear 



12 



structures significantly increases with the density at the region p > 0.16 fm~ 3 . It implies 
that more isospin-sensitive observations from heavy-ion collisions at intermediate and high 
energies and neutron stars are still required for further constraining the symmetry energy 
at the saturation and super-saturation densities. 
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